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C^l ■ Abstract 

^1^1 We are interested in the differential equations satisfied by the density of the Geometric Sta- 

■^ . ble processes G^ = {^a(*);* ^ 0}j with stabihty index a € (0,2] and asymmetry parameter 

/3 € [—1,1], both in the univariate and in the multivariate cases. We resort to their represen- 
tation as compositions of stable processes with an independent Gamma subordinator. As a 
preliminary result, we prove that the latter is governed by a differential equation expressed by 
means of the shift operator. As a consequence, we obtain the space-fractional equation satisfied 
Pm I by the density of G^. For some particular values of a and /3, we get some interesting results 

Ph . linked to well-known processes, such as the Variance Gamma process and the first passage time 

r^ ' of the Brownian motion. 
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^ ! 1 Introduction and notation 

f~^ , The Geometric Stable (hereafter GS) random variable (r.v.) is usually defined through its char- 

t:j- \ acteristic function: let Qa be a GS r.v. with stability index a G (0,2], symmetry parameter 

O I ;3 e 1-1, 1], p„s,f„n para„,eter f, £ R, scale para„,eter . > 0, then 



where 



C^ ■ (n\ _ j i — 'i/3sign{9)tan{7Ta/2), if a ^ 1 

'^"'^^^ '" I 1 + 2i^sign{e) log \6\/7r, if a = 1 

(see e.g. \9\). Moreover the following relationship holds (see [12j ) 

Eg.eef _ 1 



i-log^5£(i)W 

where 

<^^.B(e):=Ee''^^^ =exp{iefi-a''\eruj^B(e)}, eeR, (2) 

is the characteristic function of a stable r.v. Sa with the same parameters a , f3, fi, a. We will 
consider, for simplicity, the case /i = 0; then we will refer only to strictly stable r.v.'s, if a 7^ 1. 
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The main features of the GS laws are the heavy tails and the unboundedness at zero. These two 
characteristics, together with their stability properties (with respect to geometric summation) and 
domains of attraction, make them attractive in modelling financial data, as shown, for example, 
in |15j . As particular cases, when the symmetry parameter /3 is equal to 1, the support of the GS 
r.v. is limited to M"*" and its law coincides, for < a < 1, with the Mittag-Leffler distribution, 
as shown in [9] and |12) . Moreover the GS distribution is sometimes referred to as "asymmetric 
Linnik distribution" , since it can be considered a generalization of the latter (to which it reduces 
for /3 = /.i = 0, see |13j . [7|). The Linnik distribution exhibits fat tails, finite mean for 1 < a < 2 
and also finite variance only for a = 2 (when it takes the name of Laplace distribution) and is 
applied in particular to model temporal changes in stock prices (see [2]). 

The univariate GS process will be denoted as lGa{t),t > 0[ and defined by having the one- 
dimensional distribution coinciding with Qa and characteristic exponent equal to 

i>g,3 [6) = log(l + a''\9ruj^,p{9)), e G E, 

(see [21], [6]). Moreover the following representation holds 

gf(t):=5f(r(t)), t>0, (3) 

where Sa{t) is, for any t, a stable law with parameters /i = 0, /3 G [—1, 1], o" = t^' " and {V{t),t > 0} 
is an independent Gamma subordinator. We will use the following notation, for a generic process 
X :={X{t),t>0}. 

We note that, for /3 = 0, the process Qa reduces to a symmetric GS process (that we will denote 
simply as Qa), while, for /3 = 1, it is called GS subordinator (since it is increasing and Levy); we 
will denote it as Q'^. 

The space- fractional differential equation that we obtain here, as governing equations of Qa, 
are expressed in terms of Riesz and Riesz-Feller derivatives. We recall that the Riesz fractional 
derivative ^P" is defined through its Fourier transform, which reads, for a > and for an infinitely 
differentiable function u, 

T{^vy{xy,e} = -\e\''T{u{xy,e}, (4) 

where the Fourier transform is defined as J- {u{x);6} := J°°eu{x)dx (see |16j and jlOj, p. 131). 
Alternatively it can be explicitly represented as follows, for a £ (0,2], 

^Viuix) := --— -— / 1 — ^—^dz 5 

"^ ^' 2cos{a7T/2)T{l-a)dx J_^ |x - z|" ^^ 

(see [19]). The more general Riesz-Feller definition is given by 

:F{^^V^^^u{x);e} =^P1ie)TMx)■,e}, aG(0,2], (6) 

where 

^"(0) := _|^|«ei^«isne^ |^| < niin{a,2 - a} (7) 

(see |10j, p. 359 and p£]) and IP'S (6) coincides with the characteristic exponent of the stable random 

variable Sa, in the Feller parametrization, for 7 = - arctan [— /3tan ^1 . Indeed ([2]) can be rewritten 
(for jj, = 0) as 

$^^(0) =exp{cV^^(e)}, 0GM, c = a"[cos(7r7/2)]~^ (8) 

We recall now the following result on stable processes proved in [16] (in the special case c = 1), 
which will be used later: let pa{x,t), x G M, t > 0, be the transition density of the stable process 



Sa, then Pa satisfies the following space- fractional differential equation, for a E (0, 2], x G M, t > 0: 



Pa{x,0) = 6{x) 



(9) 



limM^^ p'a{x,t) = 



and the additional condition 4iPa{x,t) 



0, if a > 1. 



t=o 



Our main result concerns the space-fractional equation satisfied by the density ga{x,t), x G 
M, t > 0, of the GS process Qa- As a preliminary step we derive the partial differential equation 
satisfied by the density fr{x,t),x,t > 0, of the Gamma subordinator F and then we resort to the 
representation ([3]) of the GS process. Indeed we prove that fr{x,t) satisfies 



^ f 

ox 



-6(1 



-dt 



)/r, 



x,t > 0, 



(10) 



where b is the rate parameter of F (see (jlSp below) and e * is a particular case (for A; = 1) of the 
shift operator, defined as 



e-^^V(t) := E 



i-kdt 



n=0 



n! 



-fit) = fit -k), kGR, 



(11) 



for any analytical function / : M — )■ M. As a consequence, we show that gaix, t) satisfies, for x G M, 
t > 0, a G (0, 2], the following Cauchy problem 



^ ''^V^,pg^aix,t) = lil - e-9*)9Ux,t) 
gaix,0) = 6ix) 
lim\^\^^g^ix,t) = 



(12) 



In the n-dimensional case, we prove that the governing equation of the GS vector process in R" 
is analogous to (|12p . but the Riesz-Feller fractional derivative is substituted, in this case, by the 
fractional derivative operator V°^ defined by 



^{VXx);0} 



-i <z,e >)"M(dz) 



J'{n(x);6l}, 6I,xGM",q G (0,2], a/ 1, (13) 



where S" := {s G M" : ||s|| = 1} and M is the spectral measure (see [17) . with a change of sign due 
to the different definition of Fourier transform). The multivariate GS law has been first introduced 
in [T] (in the isotropic case) and called multivariate Linnik distribution. 

As special cases of the previous results the governing equations of some well-known processes 
are obtained: indeed, in the symmetric case and for a = 2, the GS process reduces to the Variance 
Gamma process, while, for a = 1, it coincides with a Cauchy process subordinated to a Gamma 
subordinator. On the other hand, in the positively asymmetric case, Qa reduces to a GS subordi- 
nator, which is used in particular as random time argument of the subordinated Brownian motion, 
via successive iterations (see [6], [21]) Moreover, for a = 1/2, we can obtain, as a corollary, the frac- 
tional equation satisfied by the density g'^.^ixjt) of the first-passage time of a standard Brownian 
motion B through a Gamma distributed random barrier, i.e. 

g[/2ix,t) := pjinf {S(s) = F(t)} e dx\ , x,t> 0. 



Indeed we prove that g'-^j^ix^t) the space-fractional equation 

——g[l^{x,t) = —(1 - e-^^)g'y^{x,t), x,t>0, (14) 

where d^'"^ /d\x\^'^ := ^^T^Ji , with the conditions in (fT2]) . 

2 Preliminary results 

We start by deriving the differential equation satisfied by the density of the Gamma subordinator, 
since it will be applied in the study of the equation governing the GS process (thanks to the 
representation ([3])). 

The one-dimensional distribution of the Gamma subordinator {r^ ^(t), t > 0} , of parameters 
a, 6 > is given by 

{ fe°* at-l -bx ^ > n 

^ ' ^- , t>0. (15) 

0, X < 

(see, for example, [3], p. 52). Hereafter we will consider, for the sake of simplicity, the case a = 1 
and denote Fi ^ := T. The Fourier transform of p5|) is given by 

fr{9,t):=T{Mx,ty,e}=Ee^<'^('^ = (^l-j^ , 6 eR. (16) 

Lemma 1 The density il5\) of the Gamma subordinator satisfies (for a = 1), the following equation 

^/r = -6(l-e-^*)/r, x,t>0, (17) 

where e~ * is the partial derivative version of the shift operator defined in ill]) , for k = 1. The 
initial and boundary conditions are the following 

/r(x,0) = 5(x) 
lim\^\^j^^fr{x,t) =0, t>0 



Proof. The first condition in I118\) can be checked easily by considering Iil6\) and the definition 
of the Dirac delta function, i.e. 6{x) := 2^ f^G~^^^d6. The second one is immediately satisfied by 
Iil5\). As far as equation (fJ?D is concerned, the Fourier transform of its left-hand side, with respect 
to X, is given by 



^|^/r(x,t);^l> (19) 



[bym]=-iOfT{e,t) = -id 



b -' 



b-ii 
For the right-hand side of |i7D we have that 



-bfr(,e,t) + be-9^fr{e,t) = -6 ( -A_^ + fog-^* (^_A ' 

b V / b ^'-' 



+ b 



.b-ie J \b-i9 

which coincides with <fJ9l). 



An alternative result on the differential equation satisfied by /r can be obtained by considering 
the following differential operator: for any given infinitely differentiable function f{x), 

~ (-l/kV+^ 
Ak,Jix) := Y, ■ Difix), x>0, keR. (20) 

i=i -^ 

We could use for (pOj) the formalism Ak,xf{x) = log(l + Dx/k). 

If moreover Dif{x 
follows: 



= 0, for any j > 0, the Fourier transform of (j20p can be written as 

|a;|=oo 



:F{Ak,.f{xy,d} = J2 1 / e''''Dif{x)dx (21) 

1=1 

= iog(i-|)/TO. 

Lemma 2 T/ie following differential equation is satisfied by the density of the Gamma subordinator: 

Q-^fr = -Ab,Jr, x,t>0, (22) 



with the conditions 

fr{x,0) = 6{x) 

l[mixi^^Dlfr{x,t) = 0, 1 = 0,1 



^._.. .^_n /_ni (23) 



Proof. The conditions (I23p are immediately verified by (I15p . Moreover, by taking the Fourier 
transform of the r.h.s. of ()22p . we get 



-{l«..,.} ^ l(.-f) 



-t 



= -/p(0,t)log(^l-| 

= -T{Ab,Jr{x,t);e}. 

m 

From the previous Lemma we can conclude that the infinitesimal generator of the Gamma 
process can be written as Ax = — log(l + Dx). 

3 Main Results 

3.1 Univariate GS process 

By resorting to the representation ([3]) and applying the previous results, we can obtain the differ- 
ential equation satisfied by the density of the univariate GS process Qa ■ This can be done, for 
i > 1, by considering Lemma 1 together with the result Q on 5q, as follows: by ([3]), we can write 

9i{^,t)= / P^Jx,z)friz,t)dz. (24) 



We consider hereafter the simple case 6 = 1. We then apply (|17|) . for 6=1, and we get 

(l-e-^*)5^(x,t) 

p^^ix,z){l-e-9^)friz,t)dz 





oo 







p^^{x,z){l-e-^')fr{z,t)dz 
/ P^ix,z)—friz,t)dz 

l-OO Q 

= -\p^^ix,z)fr{z,t)]Zo + J ^/^{x,z)fr{z,t)dz 

= c^^^V^J pi{x,z)fr{z,t)dz = c^^V<^^pgi{x,t). 
Jo 

In the last step we have applied the first equation in Q and we have considered that, for t > 1, 
/r(0, t) = 1. In the next theorem we prove the same result in an alternative way, which can be 
applied for any t > 0. 

Proposition 3 The density ga of the GS process Qa satisfies the following equation, for any x, t > 

and a G (0, 2], 

«^P,>^(x,t) = \l - e-^^)gi{x,t), (25) 

with conditions 

{ gi{x,Q) = 6{x) 

where c> is the spreading rate of dispersion defined in (8\). 

Proof. By (j24[) and ([8]) we can write the characteristic function of Qa as 

^^ieg^it) ^ / exp{czi;'^{e)}z'-^e-'dz (27) 



(26) 



yi-c^p'^{e)J ' 

where ipaiO) is defined in ([TD; thus the Fourier transform of the space- fractional differential equation 
(j25]) can be written as 

T{^^V^^^g^^{x,t);e} (28) 

= [by(l6D]=V^(e)^{5f(x,t);e} 



On the other hand we get 



/ \ * 

^(l-e-^*)^{gf(x,i);e} = ^(1-e-^*^' ^ > 



l-cV^°(0) 



/ \ t / \ t-i 

1 / 1 \ 1 / 1 \ 



c \ l-c^°(6i) y c I l-cV'^(6') 



which coincides with (j28|) . The conditions (|26p are clearly satisfied since 

\ t 






e 



d9 = 5{x) 
t=o 



and lim|^|^oo5a(x,t) = (by Q and ([Ml))- ■ 

3.1.1 Symmetric GS process 

In the special case of a symmetric GS process Ga we can easily derive from Proposition 3 the 
following result, which is expressed in terms of the Riesz derivative P", defined in ([3|). In its 
regularized form, for a G (0,2], the derivative ^2?" can be explicitly represented as 

^ _ r(l + g) sin(7ra/2) f^ u{x + y) - 2u{x) + u{x - y) ^ 

u^u{X) - ^ J ^^^ ay, {zy) 

(see [T6]). 



Corollary 4 The density g^ of the symmetric GS process Qa satisfies the following equation, for 
any x, t > and a G (0, 2], 

^^V^g^ix^t) = -(1 -e-^')5a(^,t), (30) 

c 

where c = a" and with conditions 

ga{x,0) = 5{x) ,^^ 

^^'^\x\^oo9a{x,t) = ■ 

Remark 5 We consider now some interesting special cases of the previous results. For a = 1, 
we show, from the previous corollary, that the density gi{x,t) of a Cauchy process C subordinated 
to an independent Gamma subordinator (i.e. the process defined as {C{T{t)),t > 0}) satisfies the 
following equation, for any x,t > 0; 



jgi{x,t) = -{l-e-''')gi{x,t), 



d\x\ 

with conditions I131\) and d/d\x\ := ^"D^.. For a = 2, we derive the governing equation of the 
density g2{x,t) of the Variance Gamma process, since the latter can be represented as a standard 
Brownian motion B subordinated to an independent Gamma subordinator, i.e. as {B{T{t)),t > 0} . 
Indeed we get that g2{x,t) satisfies, for any x, t > 0, the second order differential equation 



^ . .. 1, 
dx'^' 



^52(x,t) = -(l-e-^*)52(x,t) 
where c = a"^ and with conditions Ii31\). 



We derive now another equation satisfied by the density of the symmetric GS process, which, 
unlike (I30p . involves a standard time derivative and a space fractional differential operator which 
generalizes (|20p . Let us define the fractional version of Ak.x^ for any a > 0, as 

^IJ{^) ■■= E V ''^"'•^(^)' ^ > 0' ^ e IK' (32) 

1=1 

where ^^^ is the Riesz derivative of order i/ > 0. We note that in the non-symmetric case (i.e. for 
/3 7^ 0) we can not define the analogue to (j32]) since the Riesz-Feller derivative is not defined for a 
fractional order greater than 2. 



Proposition 6 The density g^ of the symmetric GS process Q^ satisfies the following equation, 
for any x,t > and a G (0,2], 

d 

T^ga{x,t)=A'^/^^^ga{x,t), (33) 



dt^ 



where c = a"' and with conditions 



ga{x,0) = 5{x) 

lim|xKoo^5a(2;,i) =0> / = 0, 1,... 



(34) 



Proof. The Fourier transform of (j32p is given by 



-^{^Vc,./(^);4 = Y.-^ e^'-^''vff{x)dx 



i+i 






1=1 



\og{i + c\enT{f{xy,e}. 



Therefore we get 



^{■^i/c,.9aix,t);e} = log(^^-l^)^{g„(x,t);n 



log 



i + c\e\'^J \i + c\e\'^ 

since for /? = 0, the characteristic function (|27l) reduces to 



^^iOGait) 



1 



l + cl^l^ 



(35) 



(36) 



The expression (I36p clearly coincides with the Fourier transform of the left-hand side of (j33p 

The previous result agrees with the expression of the infinitesimal generator Ax of the GS 



process, which is given by Ax = — log 



1 + 



dx^ 



(see m) 



3.1.2 GS subordinator 

In the positively asymmetric case, i.e. for /3 = 1, the process Qa reduces to a GS subordinator (we 
will denote it as Q'^). 

Corollary 7 The density g'^{x,t) of the GS subordinator Q'^ satisfies the following equation, for 
any x, t > and a G (0, 2] , 



RFtiO I 



1 



-dt\J 



'Vlig'a{x,t) = -{l-e-''^)g'^{x,t), 



where c = a'^ (cos(7ra/2)) and with conditions 



g'^{x,0) = 6{x) 



(37) 



(38) 



an 



d^^jy^i is the Riesz- Feller derivative defined by T {^^V^ iu{x);6^ = —{—i\9\)°'sign{9)T{u{x);9} 



Remark 8 We now consider the special case a = 1/2 of the previous result. It is well-known that 
the stable law with parameters a = 1/2, ^ = 0, (3 = 1, a>0 coincides with the Levy density. 
Moreover if we define as 

Tz := inf {B{s) = z} , z>0, 

s>0 

the first-passage time of a standard Brownian motion B, we have that 

P {Tz ^ dx} = p\i2{x,z), x,2;>0, 

since Tz is equal in distribution to a stable subordinator S'^j^ of index 1/2 and variance a = z^ 
(whose density is denoted as p'-^j^ix^z)). Therefore, from the previous corollary, we can derive that 
the density of the time-changed process |7r(t),i > Oj , given by 

Pl/2ix^z)fY{z,t)dz 

satisfies the following equation for any x,t > 0; 



^ jj^g[/2{^,t) = ^2(1 - e-^')g[/^{x,t), (39) 



d\x\y 



with conditions ll&8\) andd^''^/d\x\^''^ := ^^T^Ji ■ The constant in /i8y\) can be derived by considering 
that, in this case, c = -^/ct (cos(7r/4))^ and we assume that a = 1. The process Tptt) can be 
interpreted as the first-passage time of a Brownian motion through a random barrier, represented 
by a Gamma process. Thus we can conclude that 

plmf{B{s) = r{t)} edx\ , x,t>0 

satisfies the space-fractional equation i39\} . 

3.2 Multivariate GS process 

The multivariate GS distribution was first defined in p^ and applied later to model multivariate 
financial portfolios of securities, in |14] . 

In the n-dimensional case, we denote by {G"(t), t > 0} a multivariate GS process with stability 
index a G (0,2], position parameter /u = (for simplicity) and spectral measure M, then its 
characteristic function can be written as 



]ggi<e,GS(t)> 



-t 



1+1 \<9,z>\"uja,i{<e,z>)M{dz 

where S"" := {s G M" : ||s|| = 1}, < 0, z >= Y,]=i GjZj and 

1 — isign{< 6, z >) tan(7ra/2), if a 7^ 1 



(40) 



uja,i{<e,z>):- ^ i + 2isign{<e,z>)log\e\/TT, if a = 1 ' 



Moreover, as in the univariate case, the following relationship holds for the r.v. G^ := GJJ(l): 

1 
l-log$ss(^)^ 



Ee^^GS^ -^^^^, ee 



(see [14j), where 

$SsW := Ee*<^'^"> = exp{- / | < 0,z >|°w„,i(< e,z>)M{dz)}, G M" 

is the characteristic function of a stable multivariate r.v. S^ with ^ = and spectral measure M 
(see e.g. [2D], p. 65). 

Let the process {S^(t),t > 0} be defined by its characteristic function, i.e. 

$g„(^)(0) — ]Eei<e,S2(t)> ^ exp{-t / \<0,z >|"a;„,i(< e,z >)M{dz)}, 9 G M". 
Then the density p'^{x, t) of S^ satisfies the initial value problem, for a G (0, 2], a 7^ 1, 



K(x,0)=5(x) ' ^eM,t>U, 



(41) 



where c = (cos(-7rQ/2))^^ (see [T7], being careful with the signs, for the different definition of Fourier 
transform) and V^^ is the fractional derivative operator defined in ()13p . 

The results of the previous section can be generalized to the n-dimensional case, as follows. 

Proposition 9 The density g'^(x,t) of the n-dimensional GS process GJJ satisfies the following 
Cauchy problem, for a G (0, 2], a 7^ 1, 

V^,5S(x,t) = i(l-e-^05a(x,t) 

52(x,0)=,5(x) ,xGM",t>0. (42) 

lim||x|Koo52(x,i) =0 

Proof. The Fourier transform of the space-fractional differential equation in (I42p can be written 



as 



F{Vl,gl{x,t);e] 



[by ([13 



/ {-i<e,z>)"'M{dz) 



:F{g^{x,ty,e} 



< e ,z >\'^uJo.,i{< 0,z >)M{dz) 



— cos(7ra/2) 

[by (HDD] 
cos(7ra/2)(l-e-^*)-F{5^(x,t);0}, 



n -t 



1+1 \<e,Z>\''uJa,li<e,Z>)M{dz) 

gn 



by considering that 



The first condition in (j42p is verified since the characteristic function of G^, given in (j40p . reduces 
to 1 for i = 0, while for the second one we must consider that 



-i < 9,z >)" = 1 <9,z >\'^ cos{7Ta/2)uja,i{< 0,z >). 



ff2(x,i) 



p'^{x,z)friz,t)dz 



and that lim\^\_^^ p^{x, z) = 0. 
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Remark 10 If we consider the special case of an isotropic n- dimensional GS process {Gait),t > 0}, 
the previous results can be considerably simplified. Indeed in this case we can use the fractional 
Laplace operator defined by 

^{(-A)°n(x);^} = -||e||"^{u(x);0}, x,0eM" (43) 

(where \\ ■ \\ denotes the Euclidean norm) or, by the Bochner representation, as 

(_A)" = _!^^^M r^ ^"-i(^ _ A)-iAdz (44) 

(see fSf and 14J)- Moreover the n-dimensional isotropic GS process is defined through its charac- 
teristic function 

^gi<e,GJt)> ^ ( L^A Q ^n .45^ 

Vi + ||6i||"y ^ ' 

(see ii4l) CLnd its marginals coincide with the multivariate Linnik distributions introduced in |I]/. 
The process Gq can be represented as 

Gait) := S„(r(i)), t > 0, (46) 

where {Sa{t),t > 0} is an isotropic stable vector, with characteristic function 

Egi<0,s„(t)> ^ exp{-t||^||"}, e £ M". 
Then it is well-known that the density Pa{^,t) of Sa satisfies the equation 

(-A)"p„(x,t) = ^|p«(x,t) 

p«(x,0) = <5(x) , (47) 

lim||x||^ooPQ(x,t) = 

for X G M", t > 0, c = (cos(7ra/2))^^ and a E (0,2]. Therefore by Preposition 3, we can conclude 
that the density ga{x,t) of Ga satisfies the following Cauchy problem, for any x gM", t > 0; 

f (-A)"g«(x,t) = i(l-e-^*)5a(x,t) 

5c.(x,0) = (^(x) , (48) 

lim||x||->oo5a(x,t) = 
where (—A)" is the fractional Laplace operator defined in ^4^- 
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